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Cell size is a complex quantitative trait resulting from interactions between intricate genetic 
networks and environmental conditions. Here, taking advantage of previous studies that uncovered 
hundreds of genes affecting budding yeast cell size homeostasis, we performed a wide pharmaco- 
epistasis analysis using drugs mimicking cell size mutations. Simple epistasis relationship 
emerging from this approach allowed us to characterize a new cell size homeostasis pathway 
comprising the sirtuin Sir2, downstream effectors including the large ribosomal subunit (60S) and 
the transcriptional regulators Swi4 and Swi6. We showed that this Sir2/60S signaling route acts 
independently of other previously described cell size controlling pathways and may integrate the 
metabolic status of the cell through NAD^ intracellular concentration. Finally, although Sir2 and 
the 60S subunits regulate both cell size and replicative aging, we found that there is no clear causal 
relationship between these two complex traits. This study sheds light on a pathway of > 50 genes 
and illustrates how pharmaco-epistasis applied to yeast offers a potent experimental framework to 
explore complex genotype/phenotype relationships. 
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Introduction 

Cell size, similar to the majority of phenotypic characters, is 
the result of complex genetic interactions. Cell size can vary 
substantially across cell types and organisms. It is influenced 
by endogenous factors such as ploidy but also by environ- 
mental conditions. However, for a given cell type in a defined 
growth condition, cell volume distribution is constant, thus 
arguing for a homeostatic control of cell size (Jorgensen and 
Tyers, 2004). Preserving cell size homeostasis necessitates 
gauging cell size and coordinating growth (increase in 
volume) and prohferation (increase in cell number). In 
unicellular organisms such as bacteria or yeast, cell size 
drastically varies with the richness of the medium (Wright and 
Lockhart, 1965; Johnston et aU 1979), implying a supplemen- 
tary level of cell size regulation herein referred to as 'nutrient 
control'. In fact, cells grown in poor media are significantly 
smaller than those grown in rich media (Johnston et aZ, 1979) . 
Two major non-exclusive hypotheses have been proposed to 
account for nutrient control of cell size. First, decreasing cell 
size leads to an increase of the cell surface to volume ratio and 
could, therefore, be advantageous for the uptake of scarce 
nutrients (Hennaut et aU 1970; Adams and Hansche, 1974). 
Second, diminishing cell size minimizes the amount of 
biomass needed for each division round and, therefore, allows 
the increase of cell number before complete starvation 
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(Jorgensen et aU 2004; Jorgensen and Tyers, 2004). In both 
scenarios, cell size regulation by nutrients would contribute to 
an increased fitness of the population under sub-optimal 
conditions. 

In microorganisms, nutrient control of cell size has been 
studied for decades. In Bacillus subtilis (Weart et aU 2007), it 
has been shown that nutrient control of cell size occurs 
through a mechanism involving a metabolic enzyme (the 
glucosyltransferase UgtP), which, together with its substrate 
(UDP-glucose), inhibits cell division. This elegant work has 
provided the first mechanism connecting a specific metabohc 
activity to the control of cell size and prohferation. In budding 
yeast, the key notion of 'critical size' has emerged, defined as 
the minimal size required for entering a new cell division cycle 
(Hartwell et aU 1974; Johnston et aU 1977). Although the 
averaged 'critical size' is relatively constant in a defined 
culture condition, it varies with the nutrient content of the 
medium (Johnston et aU 1979). How yeast cells convert a 
'sufficient biomass signal', which could reflect volume, mass 
and/or biosynthetic capacity, into a 'division signal' is not 
entirely understood. Genetic approaches have been widely 
used to identify the key factors in this process. The first 
characterized mutants with a reduced cell size (named whi) 
affected the cyclin Cln3p. Although alleles of CLN3 stabilizing 
the protein lead to cell size diminution (Carter and Sudbery, 
1980; Sudbery etaU 1980; Nash etal 1988), knockout of CLN3 
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leads to a cell size increase (Lew et aU 1992). Cln3p interacts 
with the cyclin-dependent kinase Cdc28p and inhibits WhiSp, 
a Rb homolog that negatively regulates the MBF (Swi6/Mbpl} 
and SBF (Swi6/Swi4} transcription activators (Costanzo et aU 
2004; de Bruin et aU 2004). Thereby, Cln3p/Cdc28p stimulates 
the transcriptional activation of > 100 genes involved in the 
transition from Gl to S phase (Spellman etaU 1998). However, 
Cln3 is not essential for cell cycle progression, possibly 
because of a partial functional redundancy with Bck2 
(Ferrezuelo et aZ, 2009). The precise function of Bck2 is 
unclear, but this protein contributes to the activation of many 
genes, including most of Cln3 targets (Ferrezuelo et aZ, 2009). 
Together, these results pointed to the Gl/S transition 
machinery as a major factor in cell size regulation, yet other 
cell cycle regulators could also have a role in cell size 
homeostasis (Harvey and Kellogg, 2003). In a reciprocal way, 
cell size control contributes to Gl length variability (Goranov 
and Amon, 2010} . Moreover, cell growth capacity varies with 
cell cycle position, this capacity being higher in Gl and 
anaphase than during other cell cycle stages (Goranov et aU 
2009). Recently, Polymenis and coworkers (Hoose etaU 2012) 
further substantiated the complex relations between cell cycle 
progression and cell size control by reporting that many 
mutations disturbing cell cycle progression do not affect cell 
size. Therefore, our understanding on how the upstream 'cell 
size signals' are conveyed and integrated to the control of cell 
cycle progression remains to be clarified. 

Systematic identification of yeast cell size mutants, using 
knockout collections, has revealed the complexity of cell size 
homeostasis pathways (Jorgensen et aU 2002; Zhang et aU 
2002a} . Indeed, these authors identified hundreds of mutants 
with a median cell volume diverging significantly from that of 



the isogenic wild type. These large-scale approaches have 
revealed new master regulators (Sch9p and Sfplp) and have 
pointed to a central role for ribosome biogenesis and general 
nutrient sensing pathways (Ras and Tor) in the regulation of 
cell size homeostasis (Jorgensen et aU 2004). However, 
although important regulators have been well characterized, 
the vast majority of the identified cell size mutants, either 
small [whi) or large [Ige and nge) , have not yet been positioned 
into defined signaling pathways. Cell size control is thus a very 
interesting situation, where multiple loci contributing to a 
complex quantitative trait have been identified, but their 
organization into a network and their individual and 
combined influence remain to be elucidated. 

In this study, we used a large-scale pharmaco-epistasis 
approach, which allowed us to characterize a new pathway 
containing >50 genes and responding to a metabohc signal 
identified as NAD^ or a derivative. Effectors in the pathway 
include the sirtuin Sir2, the ribosome large subunit 60S and the 
transcription factors Swi4 and Swi6. Interaction with 
previously described master regulators, as well as intrinsic 
and extrinsic signals such as ploidy or medium richness, was 
evaluated in order to arrange this pathway within the known 
cell size controlling network. 

Results 

Cell size homeostasis is modulated by Sir2 

Searching for compounds that affect yeast median cell volume 
(hereafter referred to as 'cell size'}, we observed that 
nicotinamide (Nam) treatment of a wild-type strain resulted 
in cell size increase (Figure 1 A and B) . Nam, the amide form of 




WT sir2 sir3 sir4 bck2 cln3 Cell volume (fl) 



Figure 1 Cell size homeostasis is impaired in a sir2 mutant. (A) Wild-type (BY4742) and sir2 cell volume distributions. Strains were kept in exponential phase in 
SDcasaU medium for 48 h and then treated for 8 h with + Nam (1 00 |iM). For each strain, cell volume distributions were determined on at least 2x10"^ cells. (B) Mean 
of median cell volumes obtained for wild type (BY4742) and mutant strains grown as in A. (C) Schematic representation of Sir2 and Ysa1 enzymatic activities. 0-Ac-ADP- 
ribose and 0-Ac-ribose-5'-P stand for and 2'-0-acetyl-adenosine diphosphate-ribose and 2'-0-acetyl ribose 5'-phosphate, respectively. Protein names are written in 
red. (D) The ysa1 mutant has a wild-type cell volume. Characteristic cell volume distributions and corresponding median volumes (inset) obtained on cells grown as in A. 
Median volumes presented in B and D correspond to the mean of at least three independent determinations. Error bars indicate variations to the mean. Statistical 
analyses (B and D) correspond to an unpaired Student's f-test (***P<10"^; **P<10"^; NS, not significantly different; GraphPad Prism). 
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nicotinic acid (Na) , inhibits the sirtuin Sir2 (Bitterman et aU 
2002). Accordingly, we found that deleting S1R2 caused an 
^20% increase of the cell size, as previously reported in Yang 
et al (2011), just as did a Nam treatment on wild-type cells 
(Figure 1 A and B) . Further, Nam had no effect on 5ir2A cell size 
(Figure lA and B), thus demonstrating that Nam affects cell 
size through Sir2, most probably by inhibition of its enzymatic 
activity. As expected, the sirZ/S. large phenotype was rescued 
by the S1R2 gene reintroduced on a centromeric plasmid 
(Supplementary Figure 1} . Nam is one of the two byproducts of 
the deacetylation reaction catalyzed by Sir2, the other reaction 
product being 0-acetyl-ADP-ribose (Figure IC). The ysal 
mutant, known to accumulate 0-acetyl-ADP-ribose (Lee et aU 
2008), displayed a wild-type cell volume (Figure ID). There- 
fore, 0-acetyl-ADP-ribose does not seem to be involved in cell 
size homeostasis. 

Sir2 is a histone deacetylase involved in chromatin silencing 
(Guarente, 1999} at specific loci together— or not— with Sir3 
and Sir4. Yet, deletion of SIRS or S1R4 did not affect cell size, 
and sirS or sir4 cells were fully responsive to Nam (Figure IB), 
thus showing that the effect of Sir2 on cell size is independent 
of Sir3 or Sir4. Importantly, unlike 5ir2A, other large mutants 
such as d7i3A or bck2A were further enlarged in response to 
Nam (Figure IB), thus indicating that a maximal cell size had 
not been reached. Accordingly, Amon and coworkers 
(Goranov et al, 2009} have shown that Saccharomyces 
cerevisiae cells can reach a cell volume as big as 800 fl. We 
conclude that Nam affects cell size by phenocopying the sir2 
deletion. 



Mutants impairing ribosome biogenesis are 
epistatic to sir2 

We then wanted to determine the genetic network through 
which Sir2 acts on cell size control. In order to identify 
downstream genetic effectors, we performed a large-scale 
epistasis analysis to identify small size {whi) mutants masking 
the effect of sirZ on cell size. As the sir2 mutants is mating 
deficient (Wang et al, 2008), the combination of sir2 with whi 
mutations could hardly be done by large-scale mating and 
sporulation methods such as those developed for Synthetic 
Genetic Array (SGA) analysis (Tong et aU 2001). Instead, we 
took advantage of the fact that Nam affects cell size similarly to 
the sir2 knockout and that the effect of Nam is totally 
dependent on Sir2 (Figure lA). As Nam phenocopied the lack 
of Sir2, we used this drug to perform a pharmaco-epistasis 
analysis on 189 previously identified small size (whi) mutants 
(corresponding to the smallest mutants identified by Tyers and 
coworkers (Jorgensen et aU 2002}}. 

Among the 189 mutants, 22 mutants were clearly not 
affected by Nam (0.95 <Nam treated/untreated ratio <1.05; 
Figure 2 A red dots and Supplementary Table 1}, indicating that 
these mutants act downstream of Sir2 in the pathway. Of note, 
we found that several mutants described as whi by Tyers and 
coworkers (Jorgensen et aU 2002} were larger than wild type. 
This most probably reflects the differences in growth condi- 
tions used between the two studies. The pharmaco-epistasis 
relationship was confirmed by classical genetics. We com- 
bined sir2 deletion with 2 of the 22 mutants, rpl35b (Figure 2B} 
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Figure 2 Mutants unaffected by Nam treatment correspond to genes mainly 
implicated in ribosomal biogenesis. (A) Median volume of various whi mutants 
treated (yaxis) or not (xaxis) with Nam (100|iM) as in Figure 1A. Red dots 
correspond to Nam-unresponsive mutants (0.95 < median volume ratio + Nam/ 
- Nam < 1 .05), whereas green (WT), black (whi mutants) and blue (sch9, sfp1 
and I/I//7/5 mutants) dots correspond to Nam-responsive mutants. The 1 .23 ± 0.03 
ratio (green line) was calculated from median cell volumes obtained for three 
independent wild-type cultures. Note that the origin of both axes is set at 25 fl. 
(B and C) Mutation in either RPL35B or RAP49 genes is epistatic to sir2. 
Characteristic volume distributions were obtained on wild-type and mutant strains 
grown as in Figure 1A. Insets correspond to the mean of median cell volumes 
measured on at least four independent cultures. Error bars indicate variation to 
the mean. 
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and rpa49 (Figure 2C), and found that both mutants were fully 
epistatic to sir2 and, hence, presumably act downstream of 
Sir2 in the pathway. Another set of 74 mutants responded to 
Nam similar to wild-type cells (1.18 < Nam treated/untreated 
ratio <1.28}, indicating that the mutated gene and the Nam 
treatment probably act through independent means on cell 
size control. Interestingly, most of these mutants (42/74, 
P-value = 3.6 X 10"^^} affected various components of mito- 
chondria (Supplementary Table 2). In addition to these fully 
responsive and unresponsive mutants, the remaining 93 
mutants behaved in an intermediary way (54/93; 1.05 < Nam 
treated/untreated ratio <1.18) or appeared hyper-responsive 
to Nam (39/93; Nam treated/untreated ratio >1.28; 
Supplementary Table 2). This reveals complex gene/gene 
relationships that cannot be easily arranged into a pathway but 
may, in the future, be informative to understand the whole 
network. 

Strikingly, most of the Nam-unresponsive mutants (18/22) 
corresponded to knockout of genes encoding proteins imph- 
cated in diverse ribosome biogenesis steps (Supplementary 
Table 1). This was confirmed by GO term analyses 
(Supplementary Table 2} revealing a very significant enrich- 
ment for components of the cytosohc ribosome (13/22; 
P- value = 2.2 X 10 ~^^} and more specifically cytoplasmic large 
ribosomal subunit (10/22; P-value = 3.5 x 10 ~^^). In addition 
to the intrinsic constituents of the ribosome, six other proteins 
affected in the NAM-unresponsive mutants are involved in 
rRNA synthesis/maturation (Pihl, Uaf30 and Rpa49}, or 
required for ribosome assembly (Yvhl, Zuol and Jjjl) . Finally, 
we noticed that for two mutants, budl9 and ygrlGOw, the 
deletion of the coding region affected overlapping genes, 
namely RPL39 and NSRl, both required for ribosome 
synthesis. Hence, our pharmaco-epistasis analysis revealed a 
very strong enrichment for cytosohc ribosome mutants, 
although, quite surprisingly, strains deleted for two major 
ribosome biogenesis regulators, namely sfpl or sch9 
(Jorgensen et aU 2002; Jorgensen et aU 2004), were fully 
responsive to Nam (Figure 2 A, blue dots), indicating that Sfpl 
and Sch9 affect cell size mostly independently of Sir2. 

To get a more complete view of the role of the ribosome in 
cell size homeostasis, we measured the volume of every non- 
essential ribosomal protein mutants. Most of the yeast genes 
encoding ribosomal proteins are duphcated and, therefore, 
knockout of one of the two gene copies is generally not lethal. 
This analysis revealed a strong bias among ribosomal protein 
mutants: knockout of most of the large ribosomal subunit 
(60S} genes resulted in a Whi phenotype, whereas small 
subunit (40S} mutant cells were generally larger than wild- 
type cells (Figure 3 A and Supplementary Table 3). It should be 
stressed that for each ribosomal protein, the respective 
contribution of the two copies can be different. Therefore, 
for most of the mutants, the absence of a major cell volume 
phenotype could just reflect a minor contribution of the 
mutated gene copy. In any case, the cell volume distribution 
between 40S and 60S mutants was significantly different 
(Figure 3B} as previously observed by Hoose et al (2012), 
although in their study the significance of the difference 
between the two subunits was lower. 

The opposite effects of small- and large-subunit mutants on 
cell size may either reflect antagonistic impacts on the same 
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Figure 3 Impairment of small and large ribosomal subunits differentially affects 
cell size homeostasis. (A) Histogram of median cell volumes measured for 
mutants of the small (yellow, rpsA) or the large (blue, rp/A) ribosomal subunits 
(see Supplementary Table 3 for details). Green lines correspond to the wild-type 
mean of median cell volume (plain line) and the corresponding variation to the 
mean (dashed lines). (B) Statistical analysis of median cell volume ratios 
obtained for 40S and 60S mutants. The statistical analysis corresponds to a 
Wilcoxon test with 95% confidence intervals (***P<10"'^; GraphPad Prism). 
Wild-type median cell volumes were determined for 10 independent cultures. 
(C) DAB treatment affecting the large ribosomal subunit biogenesis leads to a 
Whi phenotype. The median cell volumes of the indicated strains grown as in 
Figure 1 and treated for 16 h with either + DAB (0.5 mg/l) or + CHX (0.01 mg/l) 
are shown. Values correspond to the mean of median cell volumes measured for 
at least three independent cultures. Error bars indicate variation to the mean. 
Statistical analysis was performed using an unpaired Student's Mest with 
95% confidence intervals (**P< 0.003; NS, not significantly different). The typical 
doubling time measured for each growth condition is shown at the bottom 
of the figure. 



pathway or separate effects on different pathways. To address 
this question, we constructed double mutants combining rpl 
and rps mutations. These double mutants showed an inter- 
mediary phenotype, suggesting that the large and small 
subunits could, at least in part, act independently on cell size 
homeostasis (Supplementary Figure 2). Together, our results 
estabhsh that the 405 and 60S ribosomal subunits do not 
contribute similarly to cell size homeostasis. Consistent with 
the genetic data, we found that diazaborine (DAB), a drug that 
specifically impairs 60S assembly (Pertschy et aU 2004), 
caused a decrease of wild-type cell size (Figure 3C}. By 
contrast, the translation inhibitor cycloheximide (CHX) had no 
major effect (Figure 3C}. Of note, in this experiment DAB and 
CHX were used at a low concentration, affecting only shghtly 
and similarly the generation time (Figure 3C). 

Together, these data point to an important role for the large 
ribosomal subunit in cell size homeostasis, most probably 
downstream of Sir2. It should be noted that those cell size 
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Figure 4 Identification of the Sir2-60S pathway downstream effectors. (A) Identification of the DAB-insensitive large (Ige) mutants. IVIedian cell volume were obtained 
on Ige mutants kept in exponential phase of growth for 36 h in SDcasaU and then treated (y axis) or not (xaxis) for 16 h with DAB (0.5 mg/l). Red dots correspond to 
unresponsive mutants (median cell volume ratio + DAB/- DAB > 0.95), whereas green (WT), black (Ige mutants) and blue (bck2, cln3, mbp1, and sir2) dots 
correspond to mutants differently responsive to DAB. The 0.83 ± 0.02 ratio (green line) was calculated from median cell volumes obtained from four independent wild- 
type cell cultures. Inset corresponds to a zoom of the region between 25 and 60 fl on both axes (purple dashed rectangle). (B) Validation of the DAB-unresponsive 
mutants by genetic epistasis relationships. Double (y axis) and single (x axis) mutants, respectively, refer to strains containing the DAB-unresponsive candidate 
mutations combined or not with the rpa49::HIS3 mutation. Median cell volumes of single and double mutants were determined on cells kept in exponential phase of 
growth for 24 h in SDcasaU medium (see Supplementary Table 5 for more details). Control (RPA49 (xaxis) and rpa49::HIS3 (yaxis)) and false-positive strains (see text) 
correspond to the green and orange dots, respectively. The 0.74 ± 0.08 ratio (green line) was calculated from median cell volumes obtained on three independent 
cultures of wild-type and rpa49 strains. (C and D) Validation by classical epistasis of the downstream role of Swi4 or Swi6 in the Sir2/60S controlling pathway. 
Characteristic volume distributions were obtained on wild-type and mutant strains grown as in Figure 1A. Insets correspond to the mean of median cell volumes 
measured on at least four independent cultures. Error bars indicate variation to the mean. 



mutants affect various steps of 60S synthesis and/or assembly, 
which take place in specific cellular compartments, suggesting 
that it is the assembled 60S particle rather than an inter- 
mediary assembled step that is important for cell size 
regulation. 



Swi4 and Swi6 act in the downstream part 
of the Sir2-60S pathway 

To further progress in our understanding of the genetic 
network involved in the Sir2/60S cell size control pathway, 
we again used a pharmaco-epistasis approach. A set of 155 
large mutants (corresponding to the largest mutants without 
major growth defect among those identified by Tyers and 
coworkers (Jorgensen et al, 2002)) were treated with DAB, 
which impairs 60S assembly and decreased cell size 
(Figure 3C). This allowed us to identify DAB-unresponsive 
mutants, potentially acting downstream of the Sir2/60S part of 
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the pathway (Figure 4A) . As the wild-type cell median volume 
ratio was 0.83 ± 0.02, we considered as unresponsive the 
40 mutants showing a ratio treated/untreated >0.95 (red dots 
in Figure 4A and Supplementary Table 4). To vahdate the 
pharmaco-epistasis approach, we used classical genetics to 
combine each of the 40 unresponsive mutants with the rpa49 
deletion, resulting in a robust Whi phenotype insensitive to 
Nam (Figure 2A) and epistatic to sir2 (Figure 2C). Among 
them, 33 were clearly epistatic to rpa49 (Figure 4B, red dots 
and Supplementary Table 5) fully validating the pharmaco- 
epistasis results, 2 were hypostatic and were considered as 
false positive (i.e., not confirmed by classical genetics; orange 
dots in Figure 4B and Supplementary Table 5) and 5 double 
mutants could not be obtained (Supplementary Table 4) . GO 
term analysis on the 33 confirmed mutants did not reveal any 
strong enrichment for a given cell component or function 
(Supplementary Table 6) that would give a clue on how the 
Sir2/60S pathway impinge on cell size. However, our attention 
was particularly drawn on swi4 and swi6, the subunits of the 
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Figure 5 Control of cell size by the Sir2/60S pathway respond to NAD^ variations. (A) Schematic representation of NAD+ synthesis and salvage pathways. Na, 
nicotinic acid (Niacin or vitamin B3); NAD + , p-NAM adenine dinucleotide; Nam, nicotinamide; NaMN, p-Na adenine mononucleotide; PM, plasma membrane. Protein 
names are in blue. (B) Na affects median volume of wild-type but not of s/r2 cells. Wild-type (BY4742) and s/r^ strains kept in exponential phase for 48 h in SDcasaU 
Na-free medium supplemented with different concentrations of Na (from 0 to 3 |iM). For each strain, the median volume measured in the absence of external Na was 
used as reference. (C and D) Intracellular NAD+ (orange bars) and Na (green bars) concentrations in wild type, s/r2 and npt1 strains. Cells were kept in exponential 
phase for 48 h in SDcasaU Na-free medium containing or not external Na (added Na, 3 |iM). Median volumes are shown in black. Metabolites were extracted and 
separated by liquid chromatography as described in the Materials and Methods section. Results presented in B-D correspond to the mean of at least three independent 
cultures for each condition and error bars indicate the variation to the mean. 



SBF complex (see Introduction section). We first verified that 
swi4 and swi6 mutants were fully epistatic to rpa49y rpl35b or 
rpl37a knockout (Figure 4C and D and Supplementary 
Figure 3 A and B). As SBF (Swi4-Swi6} and MBF (Swi6- 
Mbpl) are known to act downstream of Cln3 and Whi5 in cell 
size control, we wished to examine how DAB affected cell size 
of these mutants. By contrast to swi4 and swi6, the mbpl 
mutant showed a wild-type cell size and was highly responsive 
to DAB (Figure 4A, blue dot) . This result indicates that the SBF 
complex, but not the MBF complex, is required downstream of 
the Sir2-60S pathway. Unlike swi4 and 5u;i6, the cln3 or bckZ 
mutants deleted for the upstream effectors were responsive to 
DAB (Figure 4A, blue dot). Together, these pharmaco-epistasis 
analyses strongly suggest that Swi4 and Swi6 are positioned in 
the downstream portion of the Sir2/60S pathway, indepen- 
dently to their known upstream effectors (Cln3, Bck2 and 
WhiS). 



NAD+ as a physiological signal modulating cell 
size through the Sir2/60s pathway 

Our results estabhsh that Sir2 and the 60S subunit define a new 
pathway that contributes to cell size homeostasis. We 
subsequently questioned the nature of the physiological 
signals that would modify Sir2 activity and thereby regulate 
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cell size homeostasis. Sir2 activity is known to be modulated 
by NAD+ (Imai et al, 2000; Landry et al, 2000), which 
can be either synthesized de novo from tryptophan or 
recycled from Na (Figure 5A; Kucharczyk et aU 1998; 
Bieganowski and Brenner, 2004). Feeding wild-type 
cells with increasing concentrations of Na resulted in a 
progressive cell size decrease (Figure 5B, yellow dots). 
As expected, Na addition also caused a drastic increase in 
cellular NAD^ in both wild-type and 5ir2A cells (Figure 5C). 
However, although Na treatment affected wild-type cell 
size, this was not the case for the 5ir2A mutant (Figure 5B, 
blue dots). Importantly, a nptl mutant that is impaired for 
synthesis of NAD ^ from Na was large and did not respond to 
extracellular Na (Figure 5D, black bars). In this mutant, 
intracellular Na concentration was high and NAD^ concen- 
tration was low compared with wild type (Figure 5D), as 
expected if the recycling of Na to NAD ^ is impaired. Together, 
these results point to NAD ^ , or a derivative, as a physiological 
signal that regulates the Sir2/60S pathway. As NAD"^ is 
thought to be the natural activator of Sir2 (Imai et aU 2000) 
and as 5ir2A cells are unresponsive to Na and NAD^ 
variations (Figure 5C), we conclude that NAD^ is a very 
strong candidate as a metabohc regulator of cell size in yeast. 
We thus propose that Sir2 contributes to the control of cell size 
homeostasis in yeast in response to variations of the NAD^ 
intracellular pool. 
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Positioning of the Sir2/60S pathway in the cell size 
control network 

To get a more global view on how the Sir2/60S pathway 
impinges on cell size control in yeast, we examined the 
connections between known yeast cell size regulations and the 
Sir2/60S pathway. Cell size, similar to most complex traits, is 
not only affected by multiple genes but also by the environ- 
ment and particularly by the richness of the growth medium. 
As an example, wild-type yeast cells grown in the presence of 
raffinose, a carbon source less efficiently metabolized than 
glucose, have a 35% ±4 reduced cell volume (Figure 6A}. 
Importantly, carbon source control of cell size was still active 
in the Sir2/60S pathway mutants (Figure 6A}. Yet, the most 
downstream mutants, swi4 and swi6, which are not specific to 
the Sir2/60s pathway, are less, although still significantly, 
affected by nutrients than wild-type cells (Figure 6A}. We 
conclude that the Sir2/60S pathway acts on cell size home- 
ostasis independently of the raffinose/glucose nutritional 
control. 

In our attempt to place the Sir2/60S pathway in a more 
global network resulting in cell size homeostasis, we also 
evaluated the relationships between the Sir2/60S pathway and 
ploidy. As previously shown (Galitski et aU 1999), we found 
that cell size increased with ploidy (Figure 6B} . Yet, whatever 
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the ploidy, a Nam treatment resulted in a 27 ± 3 % increase of 
cell size, whereas DAB led to a 19 ± 4% decrease (Figure 6B}. 
From these results, we conclude that the Sir2/60S pathway and 
the ploidy control of cell size are not linked. 

Finally, we confirmed by classical genetics (Figure 6C} the 
pharmaco-epistasis relationships found for whiS (i.e., respon- 
sive to Nam; Figure 2 A, blue dot) and clnS (i.e., responsive to 
DAB; Figure 4A, blue dot) . These results imply that the Sir2/60S 
pathway modulates cell size independently of the Cln3/Whi5 
pathway, although both pathways have Swi4 and Swi6 in 
common (Figure 6D} . 

Mutations in S//72and 60S genes affect both cell 
size and replicative life span but the two 
phenomena are not strictly interdependent 

It is noteworthy that Sir2 is a major factor of the yeast 
rephcative aging, a process defined as the number of 
successive daughters produced by a cell before becoming 
senescent (Kaeberlein et aU 1999). In addition, a specific 
behavior for 60S ribosomal subunit mutants was previously 
reported in yeast for other phenotypes, including rephcative 
aging (Steffen et al 2008). The fact that SIR2, 60S mutants, 
medium richness and cell ploidy all affect both rephcative 
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aging and cell size homeostasis provocatively suggests that the 
two phenomena may be linked as previously proposed by 
others (Yang et aU 2011). To directly address this longstanding 
issue, we used combinations of mutations known to differently 
affect rephcative hfe span. As reported previously, we found a 
clear increase in the generation number for the whi mutant 
rplSla (Figure 7 A; Steffen et al 2008). However, the gcn4 
deletion, suppressing the rplSla rephcative hfe span increase 
(Figure 7 A; Steffen et al, 2008} , had no effect on the rpl31 a Whi 
phenotype (Figure 7B}. Similarly, the fohl mutation, which 
suppresses the sir2 rephcative hfe span defect (Figure 7C; 
Kaeberlein et aU 1999} , did not suppress the large phenotype of 
sir2 (Figure 7D}. Finally, as previously observed. Nam 
decreased rephcative hfe span (Bitterman et aU 2002} and this 
effect was found independent of the cell ploidy, as expected if it 
mimics a Sir2 defect (Figure 7E and F}. However, rephcative 
hfe span of untreated diploid and triploid cells was longer than 
that of haploid cells, despite the fact that increased ploidy 
resulted in an increased cell size (Figures 6C, and 7E and F} . 
Consequently, higher ploidy and Nam treatment of a wild-type 
strain both increased cell size but had opposite effects on 
rephcative life span. We thus conclude that cell size and 
rephcative hfe span, although affected by a common set of 
mutants, can be disconnected. 



Discussion 

Cell size, a highly complex trait regulated by hundreds of 
genes, offers a challenging framework to explore how complex 
genetic information is integrated in a phenotypic outcome. The 
understanding of such a complex network requires the 
identification of all the participating genes and the precise 
measurement of their individual and concerted contribution to 
the phenotype. These complex interactions between genes are 
referred to as epistasis, in the widest sense of the term. In 
yeast, epistasis is classically monitored after combining the 
mutations to be studied by mating and sporulation. This can be 
done on a large scale using approaches developed for SGAs 
(Tong et aU 2001}. However, it still requires mating and 
sporulation, which take time, and may be difficult in some 
specific cases such as for sir2 mutants, which are mating 
deficient. Here we used a pharmaco-epistasis approach based 
on the use of inhibitors to mimic specific mutations and 
thereby perform large-scale epistasis analyses. To restrict 
possible off-target effects, this approach needs the use of as 
low as possible concentrations of inhibitors, and when 
possible it should be validated by classical genetics. A major 
advantage of pharmaco-epistasis is that as no meiosis is 
required, the phenotype comparison is done on strictly 
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identical individuals in terms of genotype. Simple epistasis 
relationships are observed when one mutation (or drug 
mimicking the mutation) masks the effects of another 
mutation. We first concentrated on these situations, because 
they are relatively easy to interpret and allow positioning 
genes into linear regulatory pathways. 

In the case of cell size control, we benefit from the large set 
of data provided by Tyers and coworkers (Jorgensen et aU 
2002). We started from a set of knockout mutants showing the 
most extreme Whi (189 mutants) or Lge (155 mutants) 
phenotypes. It should be stressed that there is no theoretical 
reason to choose the most affected mutants, but this choice 
rehes mainly on practical reasons. Indeed, measuring precisely 
cell size on multiple different yeast strains is a technically 
challenging issue (Turner et al, 2012). Accordingly, in the 
course of this work, we observed that cell size is exquisitely 
sensitive to growth conditions. Even though we paid a 
particular attention to this point, the experimental error is 
routinely ± 5 % . Consequently, it is much easier and it raises 
much stronger conclusions, in terms of statistical significance, 
to study combinations of mutants with extreme phenotypes. 

We identified a subset of > 50 genes belonging to a new cell 
size homeostasis pathway (Figure 6D). Remarkably, this 
pathway includes numerous mutants affecting ribosome 
biogenesis and more specifically the large ribosomal subunit. 
Of note, it would not be surprising that sir2 affects cell size 
homeostasis through its effect on silencing of rDNA (Smith 
and Boeke, 1997). Previous work from Tyers and coworkers 
(Jorgensen et aU 2004) had estabhshed a strong connection 
between ribosome biogenesis and cell size homeostasis. Here 
we observed that small ribosomal subunit mutants tend to be 
larger than the wild-type control cells, whereas large riboso- 
mal subunit mutants are often Whi. In addition, combination 
of 40S and 60S mutants resulted in an intermediary phenotype, 
strongly suggesting that these mutants affect cell size home- 
ostasis by different mechanisms. It thus appears that there is 
not just one connection between ribosome biogenesis and cell 
size homeostasis but several layers of control. The first layer 
involves the Sfpl and Sch9 effectors on the Ribi and RP 
regulons that connect nutrient control to critical size via the 
rate of ribosome production (Jorgensen et aU 2004). Our data 
estabhsh that the Sir2/60S pathway is clearly distinct from the 
Sfpl/Sch9 network. First, as both sfpl and sch9 mutants are 
fully responsive to Nam, and second, as Sir2/60S mutants are 
responsive to nutrients, whereas sfpl or sch9 mutants are not 
(Jorgensen et aU 2004). However, the swi4 and swi6 mutants, 
the most downstream components of the pathway, although 
less responsive to the carbon source than wild-type, are not 
fully insensitive (Figure 6A). This suggests that these 
components integrate regulatory signals from separate path- 
ways (Figure 6D). A second layer of ribosomal effect on cell 
size involves the translational control of Cln3 (Polymenis and 
Schmidt, 1997). Yet, it is striking that Cln3 is neither required 
for ploidy control (Andahs et aZ, 2004) nor for nutrient control 
(Jorgensen et al, 2004) or the Sir2/60S pathway (this work). 
Further studies will be required to position Cln3 in the cell size 
homeostasis control network. The third layer, revealed by our 
pharmaco-epistasis study, specifically involves the 60S ribo- 
somal subunit. Intriguingly, a specific behavior of 60S mutants 
was previously reported in yeast for other phenotypes such as 
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ER stress response (Miyoshi et aZ, 2002; Zhao et aZ, 2003) or 
rephcative aging (Steffen et aU 2008) . However, in both cases 
the underlying molecular mechanisms remain largely 
unknown. 

In addition to the simple epistasis relationships, more 
complex gene/gene interactions were found, which cannot 
be simply interpreted. For instance, many mutants showed an 
intermediary response to Nam, indicating that they are not 
fully responsive to Nam but are neither totally irresponsive. 
Interestingly, among the 54 mutants behaving this way, more 
than half affected cytosohc ribosome. How should we interpret 
these partial effects? First, despite the fact that the mutants 
used in this study are knockout, some effects could be partial 
due to gene redundancy; this is clearly the case for ribosomal 
protein genes that are generally duphcated. Second, some gene 
interactions can be more comphcated than just phenotypic 
masking. In particular, it is likely that many of the studied 
mutants have pleiotropic effects that could simultaneously 
affect more than one pathway. In this type of analyses, one 
should keep in mind the multiple levels of complexity due to 
allehc diversity and specificity. In this perspective, yeast offers 
a simphfied model due to the possible use of knockout 
mutations and to the fact that the genetic network is analyzed 
in haploid cells. 

Another major level of complexity is interaction between 
genetic factors and the environment. Once again, yeast, as it is 
a unicellular organism, provides a simpler experimental frame. 
Strikingly, we found that the abundance of the coenzyme 
NAD^ appears as a major regulator of cell size homeostasis. 
Abundance of free NAD^ results from de novo synthesis and 
recycling from precursors such as Na, but also from its 
equilibrium with the reduced NADH form. It has been shown 
that in living cells, free NAD^ is several orders of magnitude 
more abundant than free NADH (Wilhamson et al, 1967; 
Zhang et al, 2002b). NAD^ concentration in yeast has been 
estimated to be in the millimolar range (Belenky et al, 2007; 
this work) . This concentration is by far higher than the Km of 
Sir2 for NAD + , which is about 30|iM (Bedalov et al, 2003; 
Borra et al, 2004) . Consequently, it seems likely that activation 
of Sir2 by NAD^ does not simply reflect substrate availability 
but more probably rehes on a more complex regulatory 
mechanism involving allosteric effects. Accordingly, our 
results show that moderate variations of NAD ^ concentration 
may impinge on cell size. Through its role in regulating Sir2 
activity, we propose that NAD^ acts as a 'metabohc cell sizer,' 
providing a simple and efficient mechanism to adapt cell size 
to metabohc status. 

Finally, the discovery of the Sir2/60S pathway as a major 
factor in cell size homeostasis echoed its effect on rephcative 
aging (Steffen et al, 2008), raising the hypothesis of a causal 
connection between cell size and aging. In addition, as noticed 
by Schneider and coworkers (Yang et al, 2011), many cell size 
mutants are also affected for rephcative hfe span. Accordingly, 
growth under conditions where carbon is limiting results in a 
decreased cell size and an increased rephcative hfe span. 
However, we showed here that Fobl and Gcn4, which are 
clearly involved in the rephcative aging process, do not affect 
ceh size, thus disconnecting the two phenomena. Moreover, 
both the treatment of cehs with Nam and an increased ploidy 
resulted in a larger median cell volume, while having opposite 
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effects on replicative life span. Hence, it seems that no simple 
connection between the two phenomena can yet be drawn. A 
similar conclusion was reached by Guarente and coworkers 
(Kennedy et aU 1994) on yeast cells artificially enlarged by 
a-factor arrest, although opposite results have also been 
reported (Zadrag et aU 2005). At this point, it is not possible to 
draw a definitive conclusion on this complex relationship 
between cell size and aging. Indeed, although many signaling 
pathways are shared by cell size and longevity, suggesting 
important common mechanisms, conspicuous counter exam- 
ples such as those reported in this work rather argue for a 
disconnection of the two phenomena. The fascinating ques- 
tion of relationships between size and hfe span of living 
organisms rose by Aristotle in his essay 'On Longevity and 
Shortness of Life' 23 centuries ago, remains open. 

Materials and methods 

Media, strains and plasmids 

SDcasaU is a synthetic minimal medium containing 5 % ammonium 
sulfate, 0.67% yeast nitrogen base (Difco) and 2% glucose, supple- 
mented with 0.2% casamino acids (Difco) and uracil (0.3 mM). 
SRafcasaU is a similar medium containing 2% raffinose instead of 
glucose. In specific experiments requiring modulation of Na concen- 
trations, SDcasaU was made with Na-free yeast nitrogen base 
(Formedium) and was supplemented with indicated concentrations 
of Na (Sigma-Aldrich) . The YPD medium used for replicative life span 
experiments contained 1% yeast extract, 2% peptone, 2% glucose. 
NAM (Sigma-Aldrich) and DAB (kind gift from Dr H Bergler) were used 
at indicated concentrations. All yeast strains were derived from the 
parent strains BY4741 and BY4742 of the haploid yeast ORF deletion 
collections (Winzeler a/, 1999). Polyploid strains isogenic to BY4741 
and BY4742 strains were obtained from David Pellman (Storchova 
etal, 2006). The rpa49::HIS3 mutants were obtained by transformation 
of wild-type or single mutant strains with a PGR fragment obtained on 
genomic DNA from the SL107-3B strain (Beckouet et al, 2008; generous 
gift from Michel Werner) with oligonucleotides RPA49up (5'-CGACG 
CCAATTAGCAATACTG-3') and RPA49Rv (5'-GTATTTGTAGATATGTATG 
TTGTGAG-3') . Histidine-prototrophic transformants were selected and 
insertion of the rpa49::HIS3 cassette at RPA49 locus was verified by 
PGR with oligonucleotides RPA49prom (5'-TTGTTTAGGTTGTGGGG 
TTGG-3') and RPA49Rv. All other multiple mutant strains used in this 
study were obtained by mating, sporulation and tetrad dissection. 
Double mutants were identified by PGR using the KanB oligonucleotide 
(internal to KanMX4; 5'-GTGGAGGGAGGAGGGGTAAT-3') and an 
oligonucleotide complementary to the promoter of the disrupted gene. 
The sir2 mutation leading to a decreased mating efficiency (Shore et al, 
1984), construction of double mutants was also done by crossing, but 
with a sir2 mutant strain covered by a SIR2 centromeric plasmid 
(p4099; URA3; lab collection). Meiotic sir2 segregants, having lost the 
plasmid, were then identified as uracil auxotroph. Doubling time of the 
most used strains in this study are as follows: wild-type strain 
(105 min), rplSIa (120 min), rpl35b (130 min), rpa49 (135 min), sir2 
(105 min), swi4 (165 min) and swi6 (165 min). 



Cell size distribution measurements 

All cell size measurements were performed by determining the median 
cell volume of asynchronous yeast cell cultures using a Goulter- 
counter apparatus (Beckman-Goulter) . Gompared with other methods 
used for yeast cell volume measurement, this method is fast and 
accurate, and allows to process multiple samples (Turner et al, 2012). 
Gells were grown overnight in SDcasaU or SRaffcasaU media, and 
then diluted several times in order to maintain exponential growth 
(cell number is always kept under 2 x 10 cells/ml) for 24-48 h before 
cell size measurement. It should be stressed that we noticed a 
significant effect of culture conditions on wild-type cell volume 
(Supplementary Figure 4). Gonsequently, volume comparisons 



require to be done on cells grown in the exact same conditions. To 
obtain each cell size distribution, 100 |il of culture were then diluted into 
10 ml of Isotonll (Beckman-Goulter) and size distribution of the 
population was analyzed with a multisizer4 (Beckman-Goulter) by 
counting between 10 000 and 20 000 cells for each measurement. Results 
are given as the percentage of cells counted in each of the 400-size 
classes. Median volume was obtained from the geometric cell volume 
distribution by using the Multi4 software (4.02 version; Beckman- 
Goulter) with a smoothing of 7 as previously done in yeast (Jorgensen 
etal, 2002). Statistical significance of differences between two conditions 
has been determined through the use of an unpaired Student's Mest 
(Graphpad Prism Software) . For size determination in presence of Nam 
(100 |iM), DAB (0.5 mg/1) and GHX (0.01 mg/1), cells were incubated in 
the presence of the indicated drugs for 8 h (Nam) or 16 h (DAB and GHX) 
before measurements. For epistasis studies, median cell volumes were 
determined on the four spores of at least four tetratypes. 



Replicative life span analysis 

All replicative life span experiments were carried out (at least twice) as 
described (Kaeberlein et al, 2005) on 50-100 daughters of daughter 
cells grown on standard YPD plates containing or not Nam (500 |iM, a 
dose resulting in a 20% median volume increase on plates). Statistical 
significance of replicative life span changes between strains was 
determined using a Wilcoxon rank-sum test (GraphPad Prism Soft- 
ware) using a cutoff value of P= 0.05. 



Intracellular metabolites determination 

Wild-type and mutant strains were kept in exponential phase for 48 h 
in SDcasaU Na-free medium supplemented or not with extracellular 
Na (3 |iM, corresponding to the concentration found in Na-containing 
SD medium). Metabolites extraction by rapid filtration and ethanol 
boiling method, and metabolite separation by liquid chromatography 
were performed as described previously (Hurlimann et al, 2011; 
Laporte et al, 2011). 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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